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Abstract 

We investigate the influence of the nuclear deformation on the decay rates 
of some cluster emission processes. The interaction between the daughter 
and the cluster is given by a double folding potential including quadrupole 
and hexadecupole deformed densities of both fragments. The nuclear part of 
the nucleus-nucleus interaction is density dependent and at small distances 
a repulsive core in the potential will occur. In the frame of the WKB- 
approximation the assault frequency of the cluster will depend on the geo¬ 
metric properties of the potential pocket whereas the penetrability will be 
sensitive to changes in the barrier location. The results obtained in this pa¬ 
per point out that various combinations of cluster and daughter deformations 
may account for the measured values of the decay rate. The decay rather are 
however more sensitive to the changes in the daughter deformation due to the 
large mass asymmetry of the process. 


PACS : 21.60.Gx , 23.70.+j , 27.90.+b 

I. INTRODUCTION 

The theoretical study of the heavy-cluster emission and the super-asymmetric fission 
started at the end of seventies [|IJ. Since the beginning, this phenomenon was recognized to 
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be a consequence of the shell closure of one or both fragments because of its cold nature, 
i.e. the low excitation energy involved in the process. Later on, Rose and Jones confirmed 
experimentally the existence of this new phenomenon [|2|. Since then many theoretical and 
experimental studies have been carried out (for a review see 0). Recently it was advocated 
that the cluster radioactivity is not an isolated phenomenon, and must be related to other 
processes like the cold fusion or cold fission f4j, where the closed shell effects play a dominant 
role. A still opened problem in the study of the cluster radioactivity is represented by the 
question of the existence of only the spherical or both the spherical and the deformed closed 
shells. Although both daughter and emitted cluster have in many cases, at least for even- 
even nuclei, a vanishing quadrupole deformation deformation in the ground state, higher 
multipole deformations are not railed out according to liquid drop calculations 0. Until 
now there are no experimental data available for deformed daughters. 

The influence of the ground-state deformations on the fragment emission probability have 
been considered in the past 0 . However, the inclusion of the deformation did not destroyed 
the relative satisfactory agreement already obtained between the experimental results and 
the results of simple models without deformation. The first theoretical study of the cluster 
deformation effects on the WKB penetrabilities have been carried out by Sandulescu et al. 
0 using the double folded Michigan-3 Yukawa (M3Y) potential for a spherical daughter and 
a quadrupole deformed emitted cluster. The barrier assault frequency u 0 was depending on 
the deformation only by means of the first turning point location, such that for a deformed 
configuration the barrier was shifted at a larger interfragment distance and u 0 was decreasing. 

In this paper we extend the study of the deformation effects in cluster radioactivity by 
accounting also for the deformation of the daughter nucleus and including higher multipole 
deformations, like the hexadecupole one. The interaction between the daughter nucleus and 
the cluster, in the region of small overlap and throughout the barrier is computed by means 
of a double folding potential. The nuclear part includes a repulsive core at small distances. 
In this way our deformed cluster approach supposes a cluster already formed in the potential 
pocket coming from the interplay between the Coulomb and the repulsive nuclear core on 
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one hand and the attractive nuclear force on the other hand. The depth and the width of 
this pocket will determine the assault frequency of the cluster on the barrier, through which 
it will eventually tunelate. In its turn, the penetrability will depend on the height and width 
of the barrier. Since all these geometrical characteristics depend sensitively on the shape of 
the fragments we will investigate in this paper the modification induced by the quadrupole 
and hexadecupole deformations of the fragments on the pocket and the barrier and finally 
compute decay rates for several cluster emitters. 

II. CLUSTER-DAUGHTER DOUBLE-FOLDING POTENTIAL 

The nuclear interaction between the daughter and the cluster can be calculated as the 
double folding integral of ground state one-body densities pi( 2 ){r) of heavy ions as follows 

U n (R) = j dr 1 j dr 2 Pi{r 1 )p 2 (r 2 )v(s) ( 1 ) 

where v is the NN effective interaction and the separation distance between two interacting 
nucleons is denoted by s = rq + R — r 2 , R being the centre-to-centre distance. In the past a 
G-matrix M3Y effective interaction was used to discuss light and heavy cluster radioactivity 
|7|||. This interaction contains isoscalar and isovector Yukawa functions in each spin-isospin 
( S , T) channel and an exchange component coming from the one-nucleon knock-on exchange 
term. However, this interaction is based on density-independent nucleon-nucleon forces and 
consequently very deep nucleus-nucleus potentials are obtained. A double folding potential 
based on the effective Skyrme interaction will contain a repulsive core which would prevent, 
according to the Pauli principle, a large overlap of the two interacting nuclei [9J. Thus 
the nuclear potential between two ions contains an attractive part and a repulsive one. 
Neglecting the spin dependence, it can be written as 

U n (R) = Co | Fin ((t>\ * Pi)(R) + to * pi)(R)) + itoto * AO(fl) j (2) 

where * denotes the convolution of two functions / and g , i.e. (/ * g)(x) = / f(x')g(x — 
x')dx'. The constant Cq and the dimensionless parameters Fi n ,F ex are given in ref. [9|. To 
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solve this integral we consider the inverse Fourier transform 

U n (R) = J e-^ R U N {q)dq (3) 

where the Fourier transform of the local Skyrme potential U^iq) can be casted in the form 

u N ( q ) = c 0 1 1 ' f J ' (di(<?)d2(-<7) + pi(<?)pI(-<?)) + F ex pi{q)h(q ) j (4) 

ffere p(q) and p 2 (q ) are Fourier transforms of the nucleon densities p(r) and squared nu¬ 
clear densities p 2 (r). Expanding the nucleon densities for axial-symmetric distributions in 
spherical harmonics we get 


pO) = X>A(r)Yxo(M) 


(5) 


Then 


POO 

p(q) = 4ir y i k Y xo (0 t , 0) / rdlr p.Jryj.Jqr ) 
x J » 

?(q) = 4n„(0„o) y AWfcp'*"' 2 


( 6 ) 


•o 0 0 


X 


A'A" 


r 2 drp\'{r )p\"{r )j\{qr ) 


(7) 


In this paper we take the one-body densities for both daughter and cluster as two- 
parameter Fermi distributions in the intrinsic frame for axial symmetric nuclei 


p(r) = 


Poo 


( 8 ) 


1 + exp((r — R{9))/a) 

Here poo = 0.17 fin, a denotes the diffusivity which is taken to be 0.63 for the daughter and 
0.67 for the cluster, and 


R{9) = R 0 I 1 + fo\ -^P 2 (cas9) +/3 4 ^=P 4 (cosP) 
V 47T v47T 


(9) 


is the parameterization of the nuclear shape in quadrupole /3 2 and hexadecupole (5 4 defor¬ 
mations. Here Rq = r^A 1 ^ with r 0 computed by means of a liquid drop prescription |J. 
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III. CALCULUS OF DECAY CONSTANTS 


We adopt a modified Garnmow approach || which is based on the idea that the cluster is 
pre-born, with a certain probability P 0 , in the pocket of the Skyrme+Coulomb potential and 
later on it tunnels through an essentially one-dimensional barrier. Although in the present 
case, having deformed fragments, we deal with a multi-dimensional penetration problem, we 
can disregard the orientation effects in a first approximation. It is not difficult to show that 
the value of the penetrability is maximized on the fission path corresponding to daughter 
and cluster having their symmetry axes oriented along the inter-fragment axis. For such a 
configuration the top of the barrier will attain a minimum and its location is found at a 
larger interfragment distance. 

Consequently the decay rate A will be defined as follows: 

A = uqPqP (10) 


where u 0 is the assault frequency with which the cluster bombards the walls of the potential 
pocket. It is given by the inverse of the classical period of motion 


r ° = /*y^vfyy < n > 

rtj ' 

where /j is the reduced mass of the cluster-daughter pair and r tl and r t2 are the inner turning 
points, where the potential curve intersects the Q -value (see Fig.l). Thus, in our model, z/ 0 
depends sensitively on the size of the potential pocket. The barrier penetrability is computed 
in the frame of the WKB approach 


P = exp 


dr\ 


r to. 


|(CM 


Q) 


( 12 ) 


where r t3 is the outer turning point. 

The calculus of the preformation probability P 0 is usually based on elaborated micro¬ 
scopic models. Since its calculation is beyond the purpose of this paper, we limit ourselves 
to a simple empirical formula proposed by Blendowske et.al. 0 for a, 12 C, 14 C and 16 0 
clusters 
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Po = 


(Ac < 28) 


(13) 


where the subscript c refers to the cluster and the a-spectroscopic factor is estimated as 

(P“) ewn = 6.3 x 1CT 3 and (P*) odd = 3.2 x 10“ 3 (14) 

In what follows we consider the 14 C-decay of 224 Ra. 

In Fig.l we plotted a family of potential curves U (. R ) for several quadrupole deformations 
ft? of the cluster 14 C and a fixed hexadecupolc deformation of the daughter, chosen to be 
/? 4 * = 0.008, i.e. the ground state value for 210 Pb. As one can see on this plot, the increase 
of 02 from negative to positive values is accompanied by the lowering of the barrier, while 
the bottom of the pocket goes down further. The decay rate is influenced mainly by the 
changes with deformation in the region between the first turning point and the top of the 
barrier. 

For comparison, the next plot, Fig.2, shows the variation of the interaction potential with 
the hexadecupolc deformation of the cluster. The variation brought by the hexadecupolc 
deformation is slightly different. Like in the previous case, the barrier lowers with (3 ±, but 
the bottom of the pocket rises. Obviously, this will affect the values of life-times in a different 
way than the one (3% does. The assault frequency vq will decrease with increasing cluster 
hexadecupolc deformation, but not enough to compensate the increase of penetrability. In 
overall the decay rate will increase with [3± slower than with (3% ■ 

Plots of the interaction potential U(R) for different pairs (/3-f;/^ 7 ) of quadrupole defor¬ 
mation are shown in Fig.3. One can compare the magnitude of the effect of both cluster 
and daughter quadrupole deformations and the changes which occur when we pass from 
prolate to oblate deformation. One can notice that the quadrupole deformation (3% of the 
cluster acts mainly on the right wall of the pocket, while the modification of the quadrupole 
deformation of the daughter nucleus push in opposite directions both walls of the pocket. 

The difference between positive and negative quadrupole and hexadecupolc deformations 
of the daughter nucleus can be easily understood from Fig.4. We observe different type of 
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modifications in the shape of the pocket, corresponding to fy and At deformations, respec¬ 
tively. The quadrupole and hexadecupole deformations of the daughter change the depth of 
the potential pocket in the same manner, in comparison with the case of cluster deformations 
which, as we saw in figures 1 and 2, move the bottom of the pocket in opposite directions. 

A plot of calculated A, as a function of the hexadecupole deformation of the daughter 
nucleus , for several values of (3® and a spherical cluster, is drawn in logarithmic scale in 
Fig.5. One can see, for example, that if the daughter nucleus has no quadrupole deformation 
a hexadecupole deformation /df m 0.04 may account for the experimental decay rate. Values 
of the decay rate close to the experimental one can be attained also by combinations of non¬ 
vanishing quadrupole and hexadecupole deformations. 

The next figure presents the variation of the calculated decay rate A with cluster hexade¬ 
cupole deformation for four fixed values of /3^(see Fig.6). This time the daughter nucleus 
is taken to be spherical. Recalling the observations made earlier, in connection with figures 
1 and 2, on the modification of the potential due to the quadrupole and hexadecupole de¬ 
formations one may understood why A increases faster with than with f3%. Therefore 
one may infer that it is less favorable to emit the cluster with a non-zero hexadecupole 
deformation. 

In Table I we selected some of the most favorable cases for the calculated decay rates for 
several cluster decay reactions. As we expected, prolate deformations favor the decay. In 
the case of 14 C decay of 224 Ra deformations between -0.04 and 0.04 in /Af give us decay rates 
which match the experimental one provided is taken to be not to large compared to the 
ground state value, i.e. 0.008. We may also consider a case when the cluster is deformed 
and the daughter is spherical. In the case of 14 C decay of 226 Ra, a hexadecupole deformation 
not too large than its ground state value still accounts for the experimental value although 
all other deformations are zero. 

The fact that the decay rates are more sensitive to the daughter deformations can be 
easily understood on ground of the mass asymmetry of the decaying system. An easy to 
follow explanation of this fact is given if we consider only the Coulomb part of the barrier. 
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The asymptotic part of the barrier, including the region where the third turning point r t3 is 
located, is determined essentially by the Coulomb multipoles. Introducing the quadrupolc 
moments of the charge density 


Q?’ C = fiJ 0 °°r 2 dr P ?’ c (ry 


the potential will behave for R —» oo as 

-y202\2 ZcQ2 T ZdQC 


(<?i 


000 , 


R 3 


(15) 


(16) 


Since the quadrupolc moments of the daughter and cluster depend on their 
quadrupole deformations fa it is obvious that due to the asymmetry in mass and charge 
(Zc/Zd ~ 0.1), the deformation of the daughter will have a larger influence on the barrier 
and eventually on the value of the decay rate. 


IV. CONCLUSIONS 

The aim of this paper was to extend previous studies of deformation effects in cluster 
radioactivity by considering also the deformation of the daughter nucleus and to include the 
next higher even deformation, the hexadecupole one. Considering that the cluster is pre-born 
in the potential pocket produced by the interplay between repulsive and attractive forces, 
we investigated the modifications induced by deformations on the specific potential that we 
employ. The computed decay rates depend on the assault frequency, which varies with the 
pocket depth, and on the penetrability, which changes with the barrier height. We showed 
that the experimental values can be reproduced for several selections of the deformations. If 
we maintain the cluster spherical and vary the quadrupolc and/or hexadecupole deformations 
of the daughter nucleus we may reach the experimental value within a reasonable range of 
deformations parameters. Conversely, if the daughter is spherical, then only for significant 
values of (3 we can reproduce the experimental value without taking into account . One 
might conclude that deformed states of the daughter are needed to be considered in order 



to reach the experimental A within reasonable values of cluster deformation. Configurations 
with a spherical cluster and only a hexadecupole deformation for the daughter are likely to 
occur. Also combinations of small deformations of both fragments should not be discarded. 
In conclusion, the study carried out in this paper points mainly to the importance of the 
daughter nucleus deformations, and especially its hexadecupole one. 
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TABLES 

TABLE I. Deformations of clusters and daughters which fit the calculated decay rate (see eq. (0)) to 
the experimental one \ exp . Empty spaces in the table mean null values. 


Decay 

Cluster 


Daughter 

^exp 

(s- 1 ) 

$ 


P? 

P? 

224 Ra—> 14 C+ 210 Pb 

0.160 

0.056 



9.50T0 -17 




-0.040 

0.072 






0.040 





0.040 

0.007 


226 Ra—> 14 C+ 212 Pb 




0.029 

3.50T0 -22 

234 U ^24 Ne+ 2io pb 

-0.215 

0.155 


0.098 

5.94-10' 26 


0.170 



0.010 


234 U ^28 Mg+ 206 Hg 

0.323 

-0.136 

-0.033 


2.00T0 -26 


0.113 


0.080 

0.080 


230 Th _,24 Ne+ 206 Hg 

-0.215 

0.155 

0.123 


1.60T0 -25 




0.063 
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FIGURE CAPTIONS 


Fig. 1. Variation of the interaction potential U(R) with the quadrupole deformation 
of the cluster, /3lF, when the ground hexadecupole deformation for 210 Pb is kept fixed, i.e. 
(3 4 = 0.008. The distance R is measured between the centers of mass of the two nuclei. The 
horizontal line at 30.53 MeV is the Q-value of the decay. 

Fig. 2. Variation of the interaction potential U(R) with the hexadecupole deformation 
of the cluster, fi 4 . All deformation parameters, except /3f=0.008, are zero. 

Fig. 3. Plots of the interaction potential U(R) for different pairs of daughter-cluster 
quadrupole deformation 

Fig. 4. Comparison between the effects of quadrupole and hexadecupole deformations of 
the daughter nuclei. Here, contrary to the case when the cluster is deformed, the quadrupole 
and hexadecupole deformations modify the depth of the pocket in the same manner. 

Fig. 5. The dependency of A on /3 4 in logarithmic scale for several values of daughter’s 
quadrupole deformation. The cluster is spherical in all cases. The horizontal line represents 
the experimental value for the discussed decay, \ exp = 9.50 x 10~ 17 s _1 . 

Fig. 6. Plots of the variation of calculated decay rate A with hexadecupole deformation 
of cluster for four fixed values of (3% ■ The daughter nucleus is taken spherical. 
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